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Professor H. G. Seeley’s objection to the importation of the Amphi¬ 
bian plan in explanation of a part of the Plesiosanrian shoulder girdle, 
of which other parts have been explained by reference to the Chelonian 
plan of construction, has not, I venture to think, great weight, since, 
of the early Reptilia, from the time when their remains first began 
to be studied, it has been a frequent remark that their skeletons 
comprise structural arrangements which, in existing animals, are now 
found separately. Moreover, it is thought by some of the ablest 
comparative anatomists that the Chelonian skeleton shows closer 
approach to the Amphibian than is to be found elsewhere. 


“ On Current Curves.” By Major R. L. Hippisley, R.E. 
Communicated by Major MacMahon, F.R.S* Received 
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(Abstract.) 

1. The object of the present paper is to show how to determine 
expressions for the current in circuits having iron cores, similar to 
the well-known equations 

i = — (i_ e -B(/L) 

R 


and 


i 


E 

^/(Rs+yL 2 ) 


sin (pt—6) 


for circuits without iron, which will enable the current curves to be 
pre-determined by calculation and plotted independently of experi¬ 
ment. 


d<B 

In circuits with iron cores the value of ~~~ occurring in the original 

ctt 

differential equations 


at 


(i), 


E sin pt— ( LLL =: tti... (2) 

Ctt 


continually alters as i changes. If we could obtain an expression for 
~ in terms of i, the substitution of this expression in (1) and (2) 

Cto 

should lead us to the required result. But, though such an expres¬ 
sion can be found, its substitution will generally lead to differential 
equations which cannot be solved by known methods. 

2, In the case represented by (1), where the applied E.M.E. is 
constant, we can determine by Lagrange’s formula of interpolation 
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tlie equation to the (B, H) curve of the particular core under 
consideration. This will be of the form 


B = a 0 + aiH-f a 2 W + .... +a»H w .. (3), 

where n is one less than the number of observed simultaneous values 
of B and H from which the equation is calculated ; whence 


dB 

dt 


ai L^+2 a 2 Ii 2 i- + .... +na n h<‘i“- 1 ~ 
at dt dt 


(4), 


and substituting in the equation 


we get 


E_—’ = Ri 
dt 


a Y L + 2 + 3 « 3 L 3 r -f- .... -f na n Jj >l i u ~ l 

E —Ri 


di = dt .. .. 


(5), 

(O- 


which is easily integrablo; and integrating between the proper 
limits we get after reduction 


t = l0 g Tf 


E 


E-R i 


l E E 2 E^” 1 

--aiL + ~2a 2 L 2 + —3a 3 L 3 + .. .. na }i L u 


. f l 0 E „ E 2 , E n “ 2 

i 2 a 2 L 2 +3 a 3 B° + ^ 4 aj.i +••••+ w a u Jj 


R 2 


R 3 


B n 


j . 2 / 1 j 3 E 


E 3 


E"- 


3 cl ^ ] +4 ci 4 L 4 +5 a 5 L° +••■•+ gryr 3 ncinB 1 


, / ] E E 2 E M ~ 4 

^^f^4a 4 L 4 + Si 5a 5 L 5 + s - 3 6«,L 6 + .. .. +^^1/ 


— &c., to r/y-f 1 terms ..... (7). 

The corresponding equation when the E.M. h\ is removed and the 
current is dying away is 

t z=z log — 

r 8 m 


,E 


E 2 


l E 3 


+ ^(2/ 2 L 2 - +i3/ 3 L3- 2 + j4/ 4 L 4 ^+ . 


R 


R 2 


R 3 


+ 



~~(2/,L s i+i3/ s W+i4/ 4 LV+ .... +-±-nf„IS‘i"-' S ) .. (8), 

the /i,/ 2 ,/ 3 , &c., being the coefficients of the powers of II in the 
equation to the descending (B, H) curve, which is, of course, different 
to the ascending curve. 
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3. This method is not applicable to the case in which the impressed 
E.M.F. is sinusoidal, on account of difficulties of integration. But 
both cases can be treated in another way :—Take a series of points 
on the (B, H) curve of the iron core, such that the chords joining them 
practically coincide with the curve itself. Let B K , H, c and B K + 1? H k +i 
be the coordinates of two consecutive points. The equation to the 
curve between these points is approximately 

B = m K +jll -f constant .. (9), 


where 


m K+l 


E/c + l B/C 


and therefore between these limits 


( 3/B T d'}/ / -J 

- lT = iL-_ ... (10). 

dt dt 

During the time that the current rises from i K to i K +u and B and 
H rise from B* and H* to B k+ i and H K+1 , and t rises from t K to 4-ri> 
we have 

(In 

E-m, c+1 LJ = R7 ............... (11), 


and therefore 


tx+i = t ,cH-jj— log 


E — Rf*. 
E—R^+i 


(12), 


which is true to a very close approximation for any simultaneous 
values of t and i between the above limits. From this equation, 
since t 0 and % are both zero, we can determine in succession the times 

H H 

at which the current has the known values 0, —-1,&c., 

L L 

using that value of m which applies to that particular value of H 
under consideration. In this way the current curve can be plotted. 

On making E = 0 in the original differential equation, and 
observing the proper limits, we get 


, _ , , r-r n ^L , 

hi+i — hi~\ u *■ AOg Yr 

iX £1 n+i 


( 18 ) 


as the equation to the curve representing the dying away of the 
current when the E.M.F. is withdrawn; m», m» +i being; determined 
from the descending (B, H) curve. 

Fig. 1 and Table I give the results of calculation for a circuit with 
the following constants:—Resistance, 1 ohm; E.M.F., 04315 volt; 
self-induction (without iron core), O’OOOl henry. 
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Table J. 


Ascending. ' 

Descending. 

Current. 

Time. 

Current. 

Time. 

0 *0100 

0-00023 

0 -43150 

0-000000 

0 *0200 

0 -00095 

0-26350 

0 -000235 

0-0270 

0 -00228 

0*19050 

0 -000768 

0 -0460 

0-03025 

0-14050 

0 -001460 

0*0600 

0-01649 

0-08250 

0 -003036 

0 -0825 

0-02269 

0 -04600 

0 -007670 

0-1075 

0-02696 

0 -02700 

0 -020009 

0*1400 

0-03100 

0 -01775 

0 039880 

0-1905 

0*03456 

0 -00925 

0 *067370 

0 *2635 

0 -03768 

0-00475 

0 -101070 

0-3545 | 

0 -04048 

0 -00125 i 

0*185770 

0-8750 

0 -04150 

o -ooooo 

a 

0*3875 

0 -04205 



0 -4000 

0 -04278 



0-4125 

0-04367 



0'4315 

QC 





4. When the impressed E.M.F. is sinusoidal, we substitute for 
dBjdt in the equation 

. dB . 

Esmptf*-— = Ifo .(14), 

(It 

haying determined the various values of jdt , as in the foregoing. 

As by the present method the value of m changes abruptly from 
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m K to m K+ 1 , we must employ the general solution of ( v 14), which for 
the interval 4, 4 +1 is 

* =sin +A - +ie ~ Bf/,%+iL -• (15)> 

in order that the current at the commencement of the interval 
t K , t K +i may have the same value which it had at the end of the 
interval 4_i4^ The complementary function 

A* +1 e~ Ei5 / m K+i L 

enables us to ensure this condition; for, by taking the constant 
A*+i of such a value that equation (15) is satisfied when i = i K and 
t = 4, there is no abrupt change in the current. The complementary 
function, in fact, represents the gradual dying away of whatever 
excess or defect of current there would be in the circuit when m 
changes. 

Equation (15) is true for all values of i between i K and i K +i ; and, 
therefore, enables us to find the time 4+i at which the current 
attains the known value H K + 1 /L. 

By changing ic into k + 1 we obtain similarly the time t K+2 at 
which the current has the value H K+2 /L, and so on. 

Thus the determination of 4+i is made to depend upon 4, and in 
order to make a start we must assume that the value of i is known 
for some definite value of t. When the number of alternations 
per second is not great, it is not of much consequence what assump¬ 
tion, within reason, is made, as, though the calculated curves will 
vary with the assumption made, they will all eventually merge 
into the true periodic current curve at some point which will be 
exhibited when the first evanescence of Ae~ R ^ wL takes place. 

As this complementary function is a continually decreasing quan¬ 
tity, it becomes negligible when it is allowed time enough. This 
opportunity is afforded when the straighter portions of the (B, H) 
curve are reached, and where the points on the curve can be taken 
further apart. 

When, however, the period of alternation is short in comparison 
with the time-constant of the circuit, the evanescence of Ae _E ^ mL 
' does not so readily take place; and it will generally take several 
: cycles before the current shakes down into its truly periodic forni. 
The preliminary assumption ought‘therefore in such cases to be 
made with care if it is desired to avoid the labour of calculating the 
first cycles. But, if the periodicity is quick enough, Ae~ R?/MlL can be 
taken as a constant, at any rate during the shorter intervals. 

Fig. 2 gives the plotted curve calculated for a circuit consisting of 
. 500 turns surrounding an anchor ring, having a coefficient of self- 
induction (without the core) of 0'0004i henry, and a resistance of 

s 2 
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1 ohm. The impressed E.M.F. is 12*5 volts. The periodicity is 
slow, being 6 cycles per second, and the true periodic curve is seen 
to appear before the end of the first half period. Curve A is calcu¬ 
lated with the preliminary supposition that i = 0 when t = 0, and 
Curve B on the assumption that i — 0 when pt = 0, neither of which 
is strictly correct. 

Fig. 3 gives the (B, H) curve, and Table III the observations 
upon which it is based, of the iron core under consideration. 

5. The method of this paper can be applied to the pre-determination 
of the curves of the primary and secondary currents in transformers 
and of the curve of magnetic induction with regard to time. 


[For, when there is supposed to be no magnetic leakage in the core, 
the expressions for the primary and secondary currents are— 


x 


Ey(% 4 S 2 +p 2 mVL 2 ) 

D 


sin (pt~0) -I- Ae ft 


.... (16), 


y = — ^-^ - LE eos (^_ 0)+ 


(U), 


nfpmhS 2 

w leie an — ^RS 2 -f % 2 ^ 2 m 2 L 2 (% 2 R 4- % 2 S) 9 


tan (j) — 


])mL(% 2 R 4- ni 2 S) 
~ nfES 


nf-RS 

7 “ mL(» 2 2 R+% 2 S)’ 


R being resistance of primary, S that of secondary, ni the number of 
primary turns, n 2 the number of secondary turns, L the self-induction 
of the primary. 

From (16) and (17) H, the total magnetic force on the core (being 
reduces to 

% 3 / 

H = sin (p£—0) 4 Ce~ft .... (18), 


where 


C = 


wjA4-u 2 B ^ 

% 2 


Now equation (18) treated in the same way as (15) gives the 
simultaneous values of H and t, that is, of m and t. These latter 
substituted in (16) and (17) give the simultaneous values of x and t 
and y and t. 
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If there is magnetic leakage it is necessary to discriminate between 
the H of the primary core and the H of the secondary core; these 
are respectively— 

E v '{US t +j)W(Lir-M , ) ! } sin (jpt-e) 

Ml - % ' y [ { MS -pW(LN - M 2 ) }H]> 2 mPR+LS) 2 ] 

+ C ie -y+C 2 e-y ...... (19), 


and. Eh = 


E MS sin (pt —y) 

y[{RS -p^i^LN-M-yp +pW(N R + liSp] 

+ Dje-y+:D 2 e _ y.. (20), 


1 1 

derived respectively from H a ~ — (La? + M %j) andH 2 = — (M$4Ey) ? 

»i ’ n 2 

whore * = Sln fA^+A,^ .... (21), 


and y = — 


EjpmM cos (pt—ty) 

- 


+ Btf-^4 B 2 e-v*t. ., 


.... ( 22 ). 


In these equations N is the self-induction of the secondary, M is 

the mutual induction, F is the radicle in the denominator of (19) and 

rof\\ i-i + LS 2 — p 2 m 2 N(JjN— M 2 ) 

(20), while tan « = pm —— ---—4- - -- - , 

JR(S i +jp 5i ?/i a H a ) -\-phn z M. 2 B 


tan — V m 


NR VLB 


tan g = pm 


RS-jp 2 m 2 (LN-M 2 ) 3 

L 2 S 3 -f RSM 2 -fp 2 m 2 (LlSr—M 2 ) 2 
R{LS 2 +p*m*N(LN —M 2 ) } 


ffi' = 


9> : 


• ISTRH-LS — \/{ (NR—LS) 2 4-4RSM 2 } 
2 m (LIST—M 2 ) 

NR 4 LS -f y{ (NR - LS) 2 -f 4RSM 2 } 
2m(LN—M 2 ) ! 


September } 1892,] 



Table II a .—Calculation of First Half Periods, assuming i = 0 when t 

(The * shows the point at which this coalesces with Table 115.) 
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Table lie.—Calculation of Second Half Period (i = 0 when t — 0*09184). 
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Table III.—Observations of (B, H) Curve. 1 
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* Taken from Professor Ewing’s “Magnetism in Iron, &c.,” ‘ Pliil. Trans.,’ 1885, except the first three observations, which have been 
made to show the iron starting without initial magnetism. H and B are in absolute units per sq. cm. 
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